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We report the THz response of thin fihns of the topological insulator Bi2Se3. At low frequencies, 
transport is essentially thickness independent showing the dominant contribution of the surface 
electrons. Despite their extended exposure to ambient conditions, these surfaces exhibit robust 
properties including narrow, almost thickness-independent Drude peaks, and an unprecedentedly 
large polarization rotation of linearly polarized light reflected in an applied magnetic field. This 
Kerr rotation can be as large as 65° and can be explained by a cyclotron resonance effect of the 
surface states. 



Ordered states of matter are typically categorized by 
their broken symmetries. With the ordering of spins in 
a ferromagnet or the freezing of a liquid into a solid, the 
loss of symmetry distinguishes the ordered state from 
the disordered one. In contrast, topological states are 
distinguished by specific topological properties that are 
encoded in their quantum mechanical wavefunctions [Tj. 
Frequently, a consequence of these properties is that 
there are robust "topologically protected" states on the 
sample's boundaries. The edge states of the quantum 
Hall effect (QHE) are the classic example In the 
last few years, it was realized that another class of such 
topological matter may exist in 3D band insulators with 
large spin-orbit interaction [3HS]. These so-called topo- 
logical insulators are predicted to host robust surface 
states, which exhibit a number of interesting properties 
including spin helicity, immunity to back-scattering, and 
weak aniz-localization. There are predictions of a num- 
ber of unusual phenomena associated with these surface 
states, including a proximity-efFect-induced exotic super- 
conducting state with Majorana fermions bound to a vor- 
tex [7| [S] and an axion electromagnetic response 110| , 
and proposals for applications, such as their use in tera- 
hertz (THz) devices 11 . 

Most of the signatures of topological behavior in these 
materials thus far have come from surface probes such 
as angle resolved photoemission (ARPES) and scanning 
tunneling spectroscopy [1^1 7j . These experiments have 
revealed that the surface states indeed show signatures 
of the predicted topological properties, such as a Dirac- 
like dispersion, chiral spin textures, and the absence of 
backscattering. Direct observation of the topological be- 
havior in transport has been hampered by the lack of a 
true bulk insulating state. Only recently have transport 
experiments started to distinguish the surface contribu- 
tion from the bulk [TStiST] . 

As opposed to the case of the quantum Hall effect, 
in topological insulators, the quantization of the off- 
diagonal conductivity is not a requirement for the ex- 



istence of the topological state. This, along with the 
problem of bulk conduction, has made finding a unique 
signature of this state difficult. It has been proposed 
that topological insulators may be characterized by their 
electrodynamic properties |^ due to the existence of an 
axionic term in the action AC = a0 J dxdt¥] ■ B, where 
a is the fine structure constant. is a modulo 2tt num- 
ber that distinguishes between topologically non-trivial 
{0—tt) and trivial insulators (0=0). There have been 
several proposals of how to measure this 6 parameter us- 
ing polarized electromagnetic radiation at low frequencies 
[HI • AH these proposals predict that linearly po- 

larized THz range electromagnetic waves reflected from 
(transmitted through) a sufficiently low disorder topo- 
logical insulator in a magnetic field, will undergo a Kerr 
(Faraday) rotation of the polarization plane, the magni- 
tude of which depends on the experimental configuration, 
but which is set by the fine structure constant a itself. 
Thus, such experiments can directly probe the topologi- 
cal nature of the surface states. 

In this Letter, we report the THz response of the topo- 
logical surface states (TSS) in thin films of the topological 
insulator Bi2Se3. We measured thin films of several thick- 
nesses grown on sapphire (AI2O3) substrates by molec- 
ular beam epitaxy. Due to their small thickness, low 
bulk carrier density and high mobilities, these films have 
been recently shown to exhibit thickness-independent DC 
transport [3T]. We find clear signatures of the protected 
surface states in the 2D behavior of the THz conductiv- 
ity obtained using time domain terahertz spectroscopy 
(TDTS). In addition we find a colossal Kerr rotation al- 
most entirely due to the TSS. The Kerr rotation mea- 
surement allows us to extract the effective mass of the 
2D Dirac electrons. 

In Fig. [1^ we show typical data of the real part of the 
longitudinal conductance {G^x = f^xxt) of a number of 
samples of different thicknesses {t) at 6 K and at zero 
magnetic field. The data show a clear signature of free 
electron behavior with a reasonably narrow (sa 1.2 THz 




FIG. 1. a). Real part of the conductance of three films of 
different thicknesses (16, 32 and 100 QL, 1 QL=9.4A). b). 
Spectral weight of the Drude and the phonon contributions 
obtained from the fit to the conductance (symbols). Lines 
are guides to the eye. c). Schematic of experiment, where 
multiple echoes may be separated in time (top). The bottom 
figure shows the polarizer arrangement for experiments in a 
magnetic field, the first and last polarizers can be set either 
parallel or perpendicular to each other, and the middle one 
can be rotated, d). Time domain trace of the transmission 
through a bare AI2O3 (top, displaced vertically for clarity), 
and a thin film sample (bottom) at 2 K. 



wide) Drude peak centered at zero frequency. In addi- 
tion, we observe a contribution from an optical phonon 
close to 2 THz. The data are qualitatively similar to ones 
reported by others in single crystals [271 [55] • We fit these 
data with a model consisting of two identical Drude terms 
(one for each surface), and a Drude-Lorentz term for the 
bulk phonon. These fits allow an essentially perfect pa- 
rameterization of the data (see ref.[29]) and show that the 
surface Drude term has an almost thickness-independent 
integrated spectral weight (Fig. [Tjo). The phonon spec- 
tral weight shows linear dependence with thickness, typi- 
cal of bulk response. In this fashion we conclude that the 
surface transport dominates the free electron response at 
these thicknesses. A similar conclusion has been reached 
through DC studies of these films [Hj. Using the Fermi 
velocity from ARPES [T3|, our observed scattering rate 
gives a mean free path of approximately 4.5 microns. The 
small scattering rate and long mean free path are remark- 
able considering that these films have no capping layers 
and the 'active area' has been directly exposed to atmo- 
sphere for several days. 

We now take advantage of a unique aspect of the time 
structure of TDTS and use the sapphire substrate itself as 
an optical resonator to measure both the Kerr and Fara- 
day rotation angles in the same setup (further details in 
[29j). In the experimental geometry of a thin film on a di- 
electric substrate, after transmitting through the sample. 
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FIG. 2. a). Color map of the electric field magnitude of 
the 2"'' THz pulse transmitted through the 16QL film in the 
cross polarizers configuration, for fields between -10 and 10 T 
applied along the film's c axis (positive fields are antiparallel 
to the THz pulse propagation direction), b). Real part of 
the conductance of Bi2Se3 100 QL film for magnetic fields 
between and 10 T at 5 K. 



the THz pulse partially reflects back from the substrate- 
vacuum interface and returns to the film, where it reflects 
and travels back to the detector. In principle this process 
of internal reflection inside the sapphire substrate is re- 
peated ad infinitum. These echoes are illustrated in Figs. 
[1]: and d, where we show typical time domain scans of 
a transmitted pulse and a schematic of the experimental 
configuration. In the presence of an external magnetic 
field B, the first time the pulse is transmitted the wave 
polarization may be rotated an angle ipp, the Faraday 
angle (labeled "Faraday" in fig. [1J3). As the 2nd pulse 
reflects back to the film and then reflects from it, the 
polarization may be rotated by an additional angle (pK, 
the Kerr angle (labeled "Kerr"). The fact that TDTS 
measurements are resolved in time allows the separation 
of the different contributions to the rotation angle; this 
type of separation is generally not possible with standard 
continuous wave techniques. 

In Fig. [2^ we note the first evidence of an anoma- 
lously large Kerr rotation induced by the surface states 
in applied field. The figure shows the amplitude of the 
transmitted electric field from the Kerr pulse with crossed 
polarizers PI and P3 at ±45° respectively (see [H]). As 
we increase B, we observe a rise in the transmitted am- 
plitude as the electric field apparently undergoes a large 
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rotation, since there should be no transmission for the 
crossed polarizers unless the polarization is rotated by 
the sample. In contrast, the Faraday rotation was al- 
ways small -at the level of the experimental sensitivity 
(~ 5°)- because it is proportional to the optically active 
layer thickness, which in this case is only a few nanome- 
ters. We neglect its contribution to the total Kerr -|- 
Faraday rotation of the 2nd peak in the analysis below. 
We gain insight into the origin of this apparent large Kerr 
rotation by studying the longitudinal conductance G^x in 
an applied magnetic field B. Within a conventional pic- 
ture of cyclotron resonance, we expect that the spectral 
weight of the Drude term moves to higher frequency as 
the field is increased. The shift in spectral weight is also 
accompanied by the increase in the off-diagonal conduc- 
tance Gxy with magnetic field; this increase is responsible 
for inducing rotation to polarized light impinging on the 
sample. One qualitatively observes such behavior in Fig. 
[2|d with a very weak maximum at finite frequency in the 
real part of the conductivity for the lOOQL sample at 5 K 
and for several magnetic fields up to 10 T, and with the 
rotation indicated in Fig. [2^. Films at other thicknesses 
show the same behavior. 

We quantify the rotation angle in two different exper- 
imental configurations as described in further detail in 
ref. [29) . In the coUinear polarizer configuration we ex- 
pect the amplitude of the electric field to be proportional 
to I cos(0 — if) cos{(j))\, and in the cross polarizer mode 
I cos(0 — (p) sm{(f))\, where (j> is the angle of polarizer P2 
and ip is the polarization induced by the sample. Fig. [3] 
shows precisely this behavior. Panel (a) shows a polar 
plot of the amplitude of the first and second transmit- 
ted pulses in the coUinear arrangement through a bare 
substrate at 10 T and 5 K. It is clear that no rotation 
is observed and (p=0. In Figs, and c, we show polar 
plots for both coUinear and cross polarizer geometries for 
a 16 QL sample also at 10 T and 5 K. Both sets of data 
are consistent with a rotation angle of ipK = 65±3° . This 
is an extremely large value, both in its absolute scale and 
when normalized by the field and active thickness of the 
surface state. To the best of our knowledge, it is a world 
record for the Kerr rotation of a thin film. 

In Fig. [4] we show the result of an experiment in a 
third configuration, where P2 is placed before the sam- 
ple and is rotated at a high angular speed and the in- 
and out-of-phase outputs of a lock-in amplifier give the 
electric field components X and Y (more details in j29|). 
This method allows us to perform fast scanning of the 
magnetic field and frequency dependence of the Kerr an- 
gle. In Fig. |4] we show as a function of magnetic field 
and frequency the absolute value of the Kerr angle for 
the 32 QL sample measured at 5 K. It is clear that a 
similar frequency dependence exists as in Fig. [2]3 up to 
a maximum Kerr rotation of 49.5° for this sample. The 
variation in the values of the Kerr rotation with the films 
thickness, are consistent with the variations of the carri- 
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FIG. 3. Polarizer angle (j> dependence of the transmitted elec- 
tric field amplitude normalized to its maximum value at each 
displayed frequency for the 16 QL sample. All displayed data 
were taken at 10 T and 5 K. a). Amplitude of the transmit- 
ted electric field through a bare sapphire (AI2O3) substrate 
in the coUinear polarizers configuration. Electric field of the 
Faraday pulse multiplied by 2 (open circles), and normalized 
amplitude of the Kerr pulse (closed) are shown, b). Nor- 
malized ampUtude for the Kerr pulse in the coUinear polar- 
izer configuration. A rotation angle of 65±3° is observed, c). 
Same as in b) in the cross polarizers configuration, also shows 
a rotation angle of 65±3°. 



ers concentration and mobility as found in ref. [21] 

Such a large Kerr rotation can be qualitatively ex- 
plained using the physics of cyclotron resonance, if one 
correctly takes into account the enhancement that occurs 
when the reflection off the film happens from within the 
substrate. The complex Kerr angle can be written as 

ta,n{ip ) — 2nZoGat, where n 

■n? — 1 — 2ZqGxx — Zq{GI^ + G^y) 
is the refractive index of the substrate, Zq — 377 Ohms 
is the vacuum impedance and Gxy is the Hall conduc- 
tance. We use parameters of the total conductance and 
scattering rate appropriate for our films (see |29j for de- 
tails), and obtain from the fit an effective cyclotron mass. 
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FIG. 4. Absolute value of Kerr rotation as a function 

of frequency and magnetic field for the 32 QL sample at 5 K. 
A maximum rotation of order 49.5° exists at low frequencies 
and finite fields. 



TO* ^ 0.35 nie. An estimate of the Dirac fermion cy- 
clotron mass, TO* — Ef/v'j, using the Fermi energy from 
the carrier density n ^ 3.3x lO^'^cm"^, Ef ~ 0.5eV, and 
velocity from ARPES ^14j Vf ~5xlO^ m/s, gives only 
a 10% difference from our value. We find that this for- 
mula reproduces the main features of the data in Fig. 
|4j as shown in Fig. S2 of ref. A similar analy- 

sis using bulk parameters [13 [SHI EI] cannot reproduce 
the frequency and field dependence, which further con- 
firms that the observed effect comes from the 2D surface 
states. Small values of Kerr and Faraday rotations have 
been found recently on single crystals of Bi2Se3, which 
have an appreciable bulk contribution Therefore, 
we can again conclude that the low frequency THz re- 
sponse is largely independent of the bulk contribution to 
the conductance and the observed colossal Kerr rotation 
is intrinsic to the surface metallic states. 

It has recently been found [331 [31] that another type 
of 2D state exists in crystals of topological insulators 
cleaved in ultra-high vacuum. These states originate due 
to the band-bending effects, and as a result an electron 
accumulation layer emerges at the surface; its thickness 
was found to be ~ 20 nm and the sheet carrier concentra- 
tion of the order of 10^"^ cm~^. This implies that in thin 
films of less than 40 nm, transport would appear as ef- 
fectively three dimensional. This is contrary to what we 
found in the THz conductivity, and to what is reported 
in DC transport where thickness-independent transport 
is found from 2 to 200 nm [5T] . In addition, using the ef- 
fective masses reported for these carriers (to* ~ 0.11 me 
[34]), it is not possible to reproduce the frequency and 
field dependence of the Kerr angle as shown in ref. [^ 

We have shown clear evidence for the robust THz re- 
sponse of 2D topological surface states in thin films of 



topological insulator Bi2Se3. In magnetic fields we find a 
colossal Kerr rotation with angles up to ipK ~ 65° that 
is due to the cyclotron resonance of the surface electrons. 
This measured Kerr rotation is larger by an order of mag- 
nitude than rotations typically found on high mobility 
GaAs heterostructures [32], but it is not quantized. In 
order to reach a regime where the quantized topological 
magnetoelectric effect can be distinguished from semi- 
classical cyclotron resonance physics at accessible fields, 
it appears that the films will have to have an even larger 
mobility and the chemical potential must be tuned closer 
to the Dirac point. As the Fermi energy is tuned to the 
Dirac point, measurements like ours would have the clear 
signature of a reduction of the effective Dirac mass, and 
would eventually give way to the predicted signatures 
of axion electrodynamics [S]. In general, the Kerr rota- 
tion we observe represents a benchmark for the intrinsic 
magnetoelectric effect predicted to exist in topological 
insulators. We also think that the ideas and techniques 
demonstrated in this work will be useful in studying ef- 
fects of the interaction of polarized THz radiation with 
other novel states of matter. 
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SUPPLEMENTAL MATERIAL 
Film growth. 

Thin films of topological insulator Bi2Se3 of varied 
thicknesses were grown at Rutgers by a molecular beam 
epitaxy technique on 0.5 mm thick sapphire substrates 
(AI2O3). Films were grown on ozone-cleaned surfaces 
using the two-temperature growth process. Evolution 
of the film surface during growth was monitored by 
RHEED. After deposition of 3 QL of Bi2Se3 at 110° 
C, a sharp streaky pattern was observed, indicating the 
growth of single-crystal Bi2Se3 structure. The film was 
then slowly annealed to a temperature of 220° C, which 
helped further crystallization of the film as seen by the 
brightening of the specular spot. The diffraction pattern 
and the Kikuchi lines became increasingly sharp on fur- 
ther deposition. This shows that the grown films have 
atomically flat morphology and high crystallinity. This 
process led to high quality single crystalline films with 



5 



the largest terraces, highest bulk mobilities, and lowest 
volume carrier densities as detailed elsewhere [2T| . 

THz Methods 

TDTS measurements in zero magnetic field were per- 
formed at JHU using a home-built transmission based 
time-domain THz spectrometer. In this technique, an in- 
frared femtosecond laser pulse is split into two paths and 
sequentially excite a pair of photoconductive 'Auston'- 
switch antennae on radiation damaged silicon on sap- 
phire. A broadband THz range pulse is emitted by 
one antenna, transmitted through the Bi2Se3 film, and 
measured at the other antenna. By varying the length- 
difference of the two paths, the electric field of the trans- 
mitted pulse is measured as a function of time. Ratio- 
ing the Fourier transform of the transmission through 
the Bi2Se3 film on a substrate to that of a bare refer- 
ence substrate we resolve the frequency dependent com- 
plex transmission of the film. The transmission is in- 
verted to obtain the complex conductance by the stan- 
dard formula for thin films on a substrate: T(uj) = 
[(1 -I- n)/(l + n + ZoG(a;))]e'*= where is the phase 
accumulated from the small difference in thickness be- 
tween the sample and reference substrates, n is the sub- 
strate index of refraction, Zq w 377 is the vacuum 
impedance, and G — G^^ + iG^^ is the film's effective 
complex conductance. In the case of topological insula- 
tors, the effective conductance is composed of a sum of 
a bulk contribution, Gbuik=o'buik x t where t is the film 
thickness, and a surface contribution 2Gsurface- 

Experiments in magnetic field were done in a similar 
fashion via TDTS spectroscopy at UB. The detection of 
THz radiation in this system is achieved by electro-optic 
means. In this method a ZnTe crystal is impinged upon 
by a part of the infrared femtosecond laser; when the 
THz pulse reaches the ZnTe crystal, it becomes birefrin- 
gent and changes the polarization characteristics of the 
infrared beam. The changes in the polarization of the 
infrared ultrafast laser beam, which are proportional to 
the THz electric field strength, are measured by balanced 
photodiode detection. 

As shown in Fig Ic, measurements in field were done 
in two configurations using three wire grid polarizers of 
THz radiation (PI, P2 and P3). PI was placed before the 
sample to ensure linearly polarized light was incident to 
the film, P2 and PS were placed after the sample where 
P2 was placed in a standard rotating holder that allowed 
the selection of the polarizer angle {(f>) . In the first setup 
we measured the transmitted amplitude as a function of 
P2's angle 4>, with polarizers PI and P3 parallel to each 
other at degrees (coUinear configuration). In this mode, 
we expect the amplitude of the electric field to have the 
dependence on (f) as | cos(<^ — ip) cos{(f>)\, where if is the 
pulse rotation angle measured with respect to the posi- 



tion of the first polarizer. In the second configuration, 
PI and P3 are perpendicular to each other, at 45° and - 
45° from the vertical (cross polarizer configuration) , and 
again we measure the intensity as a function of (j). In 
this case the amplitude of the electric field should have 
an angle dependence as | cos(0 — (p) sin(0)|. In the cross 
polarizer configuration, the first polarizer had an angle 
of 45° with respect to the collinear configuration due to 
the sensitivity of the electro-optic detection to the polar- 
ization of the pulse. 

The third experimental configuration used for com- 
plete characterization of the Kerr angle in the frequency- 
magnetic field plane was done using a rotating polarizer 
technique. P2 was held in a fast rotating stage and was 
placed between PI and the sample, with polarizers PI 
and P3 both collinear and oriented vertically (0°). The 
fast rotator is spun at approximately 600 rpm, and the 
resulting signal is demodulated by a lockin amplifier at 
twice the frequency set by the rotation speed. In such an 
experiment the in-phase signal of the lock-in is propor- 
tional to the transmitted electric field collinear with the 
incoming polarization (X signal), and the out-of-phase 
response is proportional to the electric field at 90° from 
the original polarization {Y signal) [35l |36] . This method 
allows fast magnetic field scanning and the analysis of the 
full spectral response. The X and Y time dependent sig- 
nals are Fourier transformed to obtain X(a;) and Y{uj). 
By taking the ratio of Y{uj) to X{uj), we obtain the tan- 
gent of the rotation angle as tan[i^(ti;)] = -|^^. In this 
work we plot the absolute value of ip. 



Temperature dependent conductance. 

We fit the zero field THz conductance data using a 
2D model dielectric constant consisting of two identi- 
cal Drude terms for the free electron response, a Drude- 
Lorentz oscillator for the phonon, and a dielectric con- 
stant (Eoo) that represents the high frequency optical 
transitions, e^x = 1 + '^^'J^^ , where a^x — Gxx/t with t 
the film thickness: 



^xx — ^oo 9 ■ "n 2 9 ■ r \^-^J 

— lujv D Up — o;^ — lujv p 

Here the subscripts D and P represent the Drude and 
the phonon contributions, respectively. We measured 3 
samples of thicknesses 16, 32, 64 and 100 QL (1 QL ^ 
0.94 nm) in the temperature range between 2 and 300 
K. Fig. |S1| shows the real and imaginary conductances 
at several temperatures for the lOOQL sample. The fit 
parameters are given in Table 1. 
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TABLE I. Fitting parameters for the Drude and phonon contributions to the conductance at zero field and 6 K. 
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FIG. SI. Frequency and temperature dependent conductance, a). Measured real and b) imaginary conductances for the 16 
QL sample for the displayed temperatures, c) and d). Respective fits using the model described in equation SI 



Calculation of Kerr angle for a thin metallic film. 

The Faraday ((/)f) and Kerr {4>k) angles are defined 



as: 



ian{ipF) ~ i 
tan((px) — i 



— r_ 
+ r_ 



where t± {r±) is the transmission (reflection) coefheient 
for (+) right circularly polarized and (-) left circularly 
polarized light. In the case of the wavelength of light be- 
ing much larger than the thickness of the film, we obtain 



the reflection and transmission coefficients as follows. 

n — 1 — ZqG± 



n 



t± = 



1 + ZoG± 
2 

1 + ZoG± 

where G± = G^x^iG^y are the (-f) right and (-) left cir- 
cularly polarized conductances. The conductances have 
the usual form for the free carrier response in a magnetic 
field as follows: 



G, 



/(47r)(r-*^) 



(r — iuj)^ — Lo'^ 

(r- 



where lOc = eB/m*c is the cyclotron frequency. 
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We can then simplify for the Kerr angle when reflection 
happens within the substrate as: 



tan((^x) = 



2nZnG 



1 — 2ZoGa: 



-^0 



(G 



2 

XX 



■ r2 ^ 

^xy J 



(S2) 



Similarly, the Kerr angle can be calculated for the case 
of reflection from the vacuum- TI interface: 



2ZoG 



xy 



l-,i2_2nZoG,,-Z2(GL + G2 ) 



(S3) 



We note that among the differences between equations 
|S2| and [S3| are the appearance of the substrate's index of 
refraction n in the numerator, which makes the rotation 
larger by this factor, as well as its multiplication in the 
denominator of the factor —2nZQGxx in eqn. |S3| which 
makes its rotation smaller (by making the denominator 
larger) . In the quantum Hall effect regime that is relevant 
for the topological magnetoelectric effect {Gxx ~ and 
Gxy ^ e^//i) these expressions can be simplified to: 



tan((p 



2nZ()Gxy 
n2 - 1 



n2 - 1 



tan((p_R-)' 



'^ZoGxy _ a 
\ — V? 1 — r 



where a is the vacuum fine structure constant. Thus, 
it is clear that making the reflection measurement from 
within the substrate leads to an enhancement of the Kerr 
angle by the substrate's index of refraction. 



With equation S2 we can reproduce the magnitude, fre- 
quency and field dependence of the measured Kerr angle 
with an almost linear dependence of the cyclotron fre- 
quency on magnetic field, as shown in Fig. S2 a), with 
an effective mass m* ^ 0.35 me obtained from the cy- 
clotron frequency in Fig. 



S2;b). 



As discussed in the 
main text, taking an effective mass of 0.12 mg, which is 
the one obtained for the states induced by band-bending 
near the surface |34) , one cannot reproduce the results of 



the experiment as shown in Fig. S2^c). Similarly poor 



comparison with the data can be obtained when the bulk 
effective mass, m* ~ 0.16 m., is used. 
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FIG. S2. Calculation of Kerr angle and cyclotron frequency, a) Calculated absolute value of the Kerr angle (Iv^kI) as a function 
of frequency and magnetic field using the parameters given in the text, b) Cyclotron frequency obtained from fitting the low 
frequency part of the spectrum (symbols), line is obtained by taking the effective mass m* ~ 0.35 me. c Calculation of the 
Kerr angle for an effective mass of m* ~ 0.12 me. 
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